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Abstract 

We study three-dimensional self-avoiding walks in presence of a one-dimensional 
excluded region. We show the appearance of a universal sub-leading exponent 
which is independent of the particular shape and symmetries of the excluded 
region. A classical argument provides the estimate: A = 2i/ — 1 ~ 0.175(1). 
The numerical simulation gives A = 0.18(2). 


1 Introduction 


An important problem in statistical mechanics is the study of the critical behaviour 
of systems in geometries with boundaries. One usually considers situations in which 
the geometric constraint changes the critical behaviour: in the renormalization-group 
terminology these are the cases in which the boundary is a relevant perturbation. 
However there are also cases, as the one we are concerned with in this paper, in 
which the leading critical behaviour is unchanged and the presence of the boundary 
appears as an irrelevant perturbation. 

Maybe because of their name the role of irrelevant operators is generally not well 
studied. Indeed, by dehnition, they do not modify the hxed point of the renormalization- 
group transformations. In particular they do not change the values of the critical 
exponents and manifest themselves only through subleading corrections to the hxed- 
point hamiltonian. Nonetheless in any actual computation the hxed-point hamilto- 
nian is replaced by an effective hamiltonian whose parameters are tuned closed to 
criticality. This means that as soon as a precise determination of the universal scal¬ 
ing behaviour is needed, it becomes important to have good control also on the terms 
responsible for corrections to scaling. 

In this paper, we will concentrate on three-dimensional self-avoiding walks (SAWs) 
in presence of an excluded one-dimensional region. We will extend the results of 

1 i i which showed for the case of a half-line the appearance of a new critical 
exponent A 0.22. Here we will consider more general one-dimensional regions and 
we will show that the value of A is independent of the shape and symmetries of the 
excluded region: only the dimensionality plays a role. Moreover we will show that A 
can be predicted to a very good accuracy by a purely geometrical argument. 

The paper is organized as follows: in Section 2 we give a general discussion of 
the relevance of the perturbation introduced by the excluded region and we give 
a prediction for the subleading exponent. In Section 3 we describe the geometries 
we analyze and the choice of observables which allow the easiest and most precise 
determination of A. In Section 4 we give a few details of the simulation, while in 
Section 5 we give the hnal results. Appendix A presents the computation of the 
generating function for ordinary random walks in presence of excluded hyperplanes 
which allow a direct check of the results of Section 2. An analogous computation 
for the case a single excluded point and for an excluded needle is reported in [Q. 
Appendix B contains some unrelated results on the small-momentum behaviour of 
the two-point function. 

2 Excluded set and corrections to scaling 

In the terminology of the held-theoretic approach to critical phenomena, the critical 
behaviour of the SAW is governed by the hxed-point of the 0{n) cr-model analytically 
continued to n = 0 i I, I, 0,1, i- The presence of an excluded region corresponds 
to a perturbation due to the introduction of an operator which creates vacancies in 
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the 0{n) model. Now consider the correlation fnnctionQ between a spin at the origin 
and one in the bnlk at location r; this fnnction will have the scaling form 

G7^(f;/3) ~ + ••• • (2-2) 


Here (3 is the inverse temperatnre, and ^ ~ (7 — is the correlation length in the 
unperturbed theory, the critical inverse temperatnre 7 and the exponent u are not 
modihed by the presence of the vacancies, nnless the exclnded region IZ is so big that 
the remaining set \ 7^ is effectively a space of lower dimensionality. See Ref. m 


for how thin a set has to be before /i = 1/7 changes. However, the behavionr of the 
other qnantities depends on whether the pertnrbation is relevant or irrelevant 0 : 


If the pertnrbation is relevant, then the leading spin-spin decay exponent rjTi 
differs from its bnlk valne rj [and as a conseqnence the leading snsceptibility 
exponent 77 ^ = (2 — rjTi)p differs from its bnlk value 7 = (2 — 7 )//]. Likewise, 
the leading scaling function Fn differs from its bulk value F; in particular, it 
has a non-trivial angular dependence . 


(b) If the perturbation is irrelevant, then 77 ^ and Fti are unchanged from their bulk 
values 7 and F. In particular, the leading scaling function F has no angular 
dependence. The effects of the perturbation show up only in the non-leading 
exponents and scaling functions 7 ^,... and F!j ^,..., which can differ from their 
bulk values. 


In either case, 77 ^ and F-ji (and indeed all of the exponents 7 ^,... and scaling functions 
Ff^,... except for an unknown amplitude) are universal in the sense that they depend 
only on the global properties of the excluded region IZ, such as its dimensionality. 

More subtle is the case in which the perturbation is marginal: rj-ji is equal to the 
bulk value but the universal scaling behaviour may be broken by logarithmic viola¬ 
tions and observables associated to the perturbation can show a complete breaking 
of universality, in the sense that they can have critical exponent with an explicit 
dependence from the coupling of the perturbation [|T^, ^ . 


To understand the effect of the introduction of the excluded region we must thus 
understand if the perturbation is relevant or irrelevant. We will resort to a geometric 
argument. Consider in a d-dimensional space a set E and let d{E) denote the number 
of dimensions in which the set E extends to inhnity. Let us now recall the fundamental 
rule in geometric probability for the dimension of the generic intersection H ft R of 
two geometric sets A and B which are immersed in a space of d dimensions: 


d{AnB) = d{A) + d{B)-d . 


(2.3) 


^In the SAW language we have 


( 2 . 1 ) 

N=0 


where CN,Tz{r) is the number of SAWs going from 0 to r in A^ steps without intersecting the region 
7 ^. 
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A negative sign of d{A fl B) means that A and B do not generically intersect in 
dimension d outside any bounded volume. The application of (|2.3|) extends also to 
random geometries where one considers a probability measure on a conhguration 
space which is concentrated on a set of events with given Hausdorff dimension. 

For example, the generic intersection of two ordinary random walks, which have 
Hausdorff dimension 2, has dimension 4 —d, and thus they do not generically intersect 
in d > 4. By using the well-known random-walk representation of the euclidean 
0 (n)-vector held theory [^, H, for which the interaction is concentrated on the 


intersections among walks, it is then possible to understand why for d > 4 in the 
critical region only a trivial theory is recovered: simply because the walks intersect 
almost nowhere! We can say that for d > 4: 

a) the probability of intersections of two random walks in the critical region scales 
towards its limiting value with the correlation-length as 

b) critical indices take their free-held (that is mean-held) values, the interaction 
term in the hamiltonian is irrelevant and induces only sub-leading corrections. 

In d = 4 the interaction becomes marginal and is responsible only for logarithmic 
corrections to the critical indices 0]. It is said that 4 is the upper critical dimension 
of the model. 

In d < 4 the interaction is relevant and thus it changes the critical indices. 

Similar ideas have been used to discuss the critical behaviour of gauge theories 
and random surfaces, see for example []^ . 


Let us now consider the case of walks in presence of an excluded region TZ of 
dimension d-ji. Consider hrst ordinary random walks whose Hausdorh dimension is 
two. Then the dimension of a generic intersection with the region IZ of dimension d^ 

is 

dint = 2 -|- d-ji — d (2.4) 


The previous argument suggests the following cases 


• dji = 0: this is the case in which we are excluding only a hnite set of lattice 
sites. The upper critical dimension (with its logarithmic corrections) is d = 2; 

• d-ji = 1: this is the case in which we are excluding a hnite set of one-dimensional 
lines. The upper critical dimension (with its logarithmic corrections) is d = 3. 

Formula (p.4|) can also be obtained from other considerations [^j. If 7?. is a dji- 
dimensional hyperplane, consider the projection of the walk on TZ', the orthogonal 
complement of TZ. Also the projection is a random walk. The probability of in¬ 
tersection will be then the probability of hrst return to the point TZ fl TZ'. But the 
probability that a d-dimensional random walk eventually passes through a point scales 
as l/^d/ 2 -l if ^ 

> 2, while it is one if d < 2. Since in our case d = d — dti, we see 
that for dint > 0 the walk generically intersects the region TZ, while for dint < 0 the 
probability vanishes as agreement with our argument. Finally notice 

that when dji = d — 1 , one can apply the results for the so-called vicious walkers. 
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Indeed the behaviour of d one-dimensional vicious walker^ can be viewed as the be¬ 
haviour of a d-dimensional walker in a space bounded by d{d—l)/2 (d—l)-dimensional 
hyperplanes |^, 


Let us now arrive at the case of SAW. Their Hausdorff dimension is 1/n (see for 
example 0, Appendix B]). For the dimension of a generic intersection with the region 
TZ we get 

dint = — I" dti — d (2.5) 


V 

Let us consider once more a series of cases, remembering that for d < 4 according to 
the Flory formula z/ 3/(d -|- 2), which is exact in d = 1,2, must be corrected by 
logarithmic violations in d = 4, and is a good approximation in d = 3. 


• d 7 ^ = 0: Excluded points are a relevant perturbation only if d < 1/z/, that is 
in d = 1. We remark that the sub-case in which the excluded set consists of a 
single point TZ = V, when V is chosen to be a nearest-neighbour of the origin 
of the walk, can be mapped into a problem without vacancies. Indeed each 
walk of N steps, starting from the origin, can be seen as a walk of A^ -|- 1 steps 
starting from V whose hrst step is the previous origin. Under this mapping 
the asymmetry induced by the exclusion of V can be seen as the correlation 
between the position of the end-point of the walk with the direction of the hrst 
step. 

• d-ji = 1: in d = 1,2 the perturbation is relevant, but it is already irrelevant in 
d = 3. 


Thus in d = 3 with a hnite set IZ of excluded lines, we expect that, regardless of 
their disposition is space, the probability that a walk intersects the region TZ scales 
as , or, since ^ ~ A^^, as . 


In the held-theoretic language this dimensional argument implies in (p.2|) 


Tj' = 7]- dint = y] + 2- 




( 2 , 6 ) 


Consequently, if Pk{x) is a homogeneous polynomial of degree k in the coordinates of 
the end-point of the walk, we have 


{Pk{x))N,n = iV"" 


A{k) + 


Bn{k) 

ATA 



(2.7) 


where ( • ) is the average in the ensemble of walks of length N that do not intersect 
IZ and A a sub leading exponent. If (|2.6|) holds, we have 

A = -dintk' = 2 z/ - 1 . (2.8) 

In general we expect renormalization effects to change (p. 6 |) and thus (p. 8 |) intro¬ 
ducing an anomalous dimension. However experience with three-dimensional models 

^ Related work on vicious walkers can be found in Ref. ^ . 
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indicates that these corrections, if not vanishing, are extremely small and thus we 
expect ( p.8|) to be in any case a very good approximation. 

The amplitude A{k) is not universal, because of an unknown scale factor, but 
it should not depend on the excluded region. Moreover, since in the critical limit 
the distribution of the end-point is rotationally symmetric, A{k) vanishes whenever 
/ dQxPk{x) = 0 where dflx is the normalized measure on the sphere This 

suggests a very convenient way to compute the subleading exponent A induced by 
the introduction of the excluded region. The idea is to consider observables which 
have zero expectation value in the rotationally-invariant continuum limit and which 
do not vanish under the residual discrete symmetry that the lattice has after the 
introduction of the excluded region. For these quantities A{k) = 0 so that the leading 
term scales as which makes the determination of A much easier. Let us notice 

that in all this discussion we have always assumed that the hrst subleading exponent 
is related to the introduction of the excluded region. This is true for those quantities 
for which {Pk{x)) would vanish in absence of the excluded region, i.e. for those Pk{x) 
which are not only non-rotationally invariant but which are also not symmetric under 
the transformations of the cubic group. If instead {Pk{x)) would not vanish even in 
absence of the excluded region, it is not clear a priori which exponent should show 
up hrst. For cubic-symmetric Pk{x) for which A{k) = 0 an extensive analysis [^, 
shows that in absence of any excluded region 


{Pk{x)) ~ A 


ku-Alatt 


(2.9) 


with Aiatt ~ 2z/. For d = 3 this exponent is much larger than (|2.8|) and thus, in 
presence of an excluded region, we expect all non-rotationally invariant Pk{x) to 
behave as with A given by (|2.8|) . We should also mention that the exponent 

A dehned here is unrelated to the exponent uju ~ 0.5 which controls the leading 
conhuent correction for SAWs in absence of any excluded region and which has been 


the object of much attention in the literature ^ ^ . 


3 The models 

We have concentrated upon the problem of 3-d SAWs starting from the origin in 
presence of an excluded region TZ which consists of a hnite collection of half-lines — 
we will call them ''^needled' — along the coordinate axes. 

For this purpose we have studied the following cases (see hgure P for the region 
TZ (for a reason to be clarihed later all the needles start from a site whose distance 
from the origin is two: for example the needle along the positive 2 ;-axis will start from 
the site with coordinate (0,0,2)): 

1 ) a needle along one direction; 

2 ) two needles along the positive and negative directions of an axis; 

3) two needles along two different axes; 
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4) four needles in a plane along the two directions of two axes; 

5) three needles along three different axes; 

6 ) six needles along the two directions of all the three axes. 

Notice that in the last case it is necessary that the needles start from a point 
located at distance two from the origin, otherwise the SAW starting from the origin 
would necessarily touch one of the needles. 

Let us now dehne a few observables which will allow us to study the effect of 
the introduction of the excluded region. If (x, y, z) are the coordinates of the end¬ 
point of the SAW (or, equivalently, (r, 0, 0) in polar coordinates) a natural choice of 
observables is given by 0)) ^ 7 ^ 0, where 0) are the spherical harmonics. 

If (rW; „j(^,0)) is not identically vanishing because of the residual cubic symmetry 
which survives the introduction of the excluded region, this quantity is a natural 
candidate for a direct determination of the exponent A. Indeed, since the critical 
limit is rotational-invariant, the leading term for N —> cx), i.e. A{k) in formula (|2.7|) , 
will vanish and thus this quantity will scale as 

Let us now classify the various observables according to the value of 1. For I = 1 
the possibilities are 


r(W,i(0,0) + W,_i(0,0)) ~ X (3.1) 

r(Fi,i(0,0)-W,_i(0,0)) ~ y (3.2) 

rYi,o(6',0) ~ (3.3) 

We can thus dehne Oi{x) = x and the obvious permutations. This observable is 
useful whenever the region TZ is not invariant under the inversion of the x-axis. 

Let us now consider 1 = 2. In this case we have 

r‘^{y2,2{d,(t))+ Y2-2{0,(f))) ~ x'^-y'^ (3.4) 

r\Y2,2{e,4>)-Y2,-2{^.4>)) ~ xy (3.5) 

r\Y2,i{0,4>) + Y2,-i{0,4>)) ~ (3.6) 

r'^{Y2,i{d,4>)-Y2 -i{9,4>)) ~ yz (3.7) 

r%o(0,0) ~ x^ + y^-2z^ (3.8) 

We can thus dehne three observables 

02,i{x,y) = x^-y^ (3.9) 

02,2{x,y) = xy (3.10) 

02,3{x,y) = (3.11) 


It is clear that, whenever a symmetry between two axes exists, the two variables 02,i 
and 02,3 are equivalent. In our calculation we have thus only considered the last two 
quantities. 
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Geometry 1 


Geometry 3 



Geometry 5 


Geometry 2 




Geometry 4 



Geometry 6 


Figure 1: The six different geometries we have simulated. For each geometry we 
have reported as a dot the starting point of the walk (located at the origin) and the 
starting points of the excluded needles (located at distance two from the origin). 









We have finally considered / = 4. In this case many different observables can be 
deffned. We have only considered 

04,1 = ^ x'^+ y'^ + z'^+ x^z^+ y‘^z^) (3.12) 

Oa, 2 {x) = X^ + ]^{y^ + - 3x^ (y^ + z^'j (3.13) 

Notice that when the excluded region is symmetric under any exchange of two axes 
O 41 and O 42 are equivalent. It is also worth remarking that O4 1 and O42 do not 
vanish even in absence of any excluded region, as they are invariant under the full 
cubic group. 

Let us now define exactly our observables for the various geometries: 

1. geometry 1 [needle along the positive x-axis]: we measure Oi = Oi{x), O 2 = 
02 ^z{y,z) and O4 = 04,2(2;); 

2. geometry 2 [two needles along the x-axis]: we measure O2 = 02 , 3 {y,z) and 
O4 = 04,2(2;); 

3. geometry 3 [two needles along the positive x- and |/-axis]: we measure Oi = 
l{ 0 i{x)+ 0 i{y)), O 2 = ^{02,3{y,z)+02,3{z,x)) = -^ 02 ,six,y), O2 = 02,2(2;,?/) 
and O4 = 04,1; 

4. geometry 4 [four needles in the 2;?/-plane]: we measure O2 = ^{02^3{y,z) + 
02,3(2:,2;)) = -l 02 , 3 ix,y) and O4 = ^(04,2(2;) + 04,2(1/)); 

5. geometry 5 [three needles along three different axes]: we measure Oi = |(0i(2;) + 
Oi{y) + 01(2:)), O2 = 1(02,2(2;,?/) + 02,2(2;, 2:) + 02,2(2;,?/)) and O4 = 04,1; 

6 . geometry 6 [six needles]: we measure O4 = 04,1. 

4 The Monte Carlo simulation 

4.1 Monte Carlo observables 

The purpose of our simulation was to compute the mean values of the observables we 
have deffned in the previous Section in presence of an excluded region TZ. A direct 
strategy would be to simulate SAWs in presence of the region IZ and then to compute 
the mean values of the various observables in the usual way. However this strategy 
requires different simulations for different excluded regions TZ. To avoid repeating the 
runs many times we have simulated SAWs without any excluded region and then we 
have reweighted the results in order to obtain the mean values of interest. Since, as 
we have previously discussed, the introduction of the excluded region is an irrelevant 
perturbation, this strategy does not introduce any significant loss of efficiency. 
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Let us suppose we want to compute {0)n,ti where {■)n,'R indicates the ensemble 
of SAWs of length N that do not intersect the region TZ. Then we simply use 






(4.1) 


where {■)]sf indicates the average in the ensemble of all SAWs of length N and 6 ^ is 
an observable which assnmes the value one if a walk does not intersect TZ and zero 
otherwise. Notice moreover that gives also the probability pn that a SAW 

intersects the excluded region as = 1 — {9^'^'^)^. 


We have used (^T|) to obtain from our simulation the mean values indeed 

it is enough, beside measuring at each Monte Carlo step i the value Oi, to record 

('R.) 

91 which assumes the value one if the walk intersects TZ, zero otherwise and then 
estimate {0)]\f^Ti by 






(4.2) 


l^i=l 

In order to further reduce the variance of our estimates we have also symmetrized 
our observables. To clarify the method suppose we want to compute {Oi)N,n for 
geometry 1 (JZ is the positive x-axis). We could consider 


■ 

A symmetrized alternative is 

+ y[e^+y) - e^-y^) + z{9^+^^ - 9 ^-^^))n 

(^Qi+X) ^ Q{-x) Q{+y) Q{-y) Q{+z) ^ 6l(-d)^ 


(4.3) 


(4.4) 


It is obvious that this quantity has the same mean value of the previous one. However 
this symmetrization reduces the error on the estimates essentially at no computational 
cost. Indeed, since we study all six geometries at the same time, we must check in 
any case if the walk intersects any of the axes. 

In general, given a geometry TZ we have considered all the possible sets sti of 
excluded needles which can be obtained by all the transformations of the cubic group 
(six for geometry 1, three for geometry 2, twelve for geometry 3, three for geometry 4, 
eight for geometry 5, one for geometry 6). To each region sti we associate a variable 
which is one if the walk does not intersect the set st^, zero otherwise and 
which is the suitably transformed variable. Then we compute from 


{0)N,n = ^ 


(4.5) 


The symmetrization has a large effect on the static variances of the various observ¬ 
ables. For instance consider for each geometry the ratio 


R 


var 


var 9^ 

J- y 0i^n) 

n-R. ^ 


(4.6) 
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where “var” indicates the static variance and n-ji is the number of terms in the sum. 
We hnd R = 6.6, 3.2, 3.8, 2.7, 1.8, 1 in the six geometries, which represents a 
considerable improvement. The improvement on the error bars is however not so 
large as the symmetrized observables are more correlated than the unsymmetrized 
ones. Thus the error bars are only 20-30% better. 

To conclude let us comment briefly on the determination of the error bars. From 
a Monte Carlo run of n iterations we have estimated (O) using 


{0)n,71 — 


q{s-r.)q{stz} 

Wl 


(4.7) 


where, as usual O = - To compute the error bars one must take into account 

the correlation between the numerator and the denominator: the independent error 
formula indeed overestimates the true error bar. For Oi the difference is of about 
20—30%, for O 2 and O 2 of 4—7%, while for O4 the difference is negligible. We have 
thus used the following relation, valid in the large sample limit. 


Var 


where A is given by 


( q{s-r.)q{s-r)\ 

Wl 

V {^7j} y 


(E*' 

{■5^} 


s-R.) Qi^Tz)^ 


{Y. 

{«7^} 


-Var (M) 


(4,8) 


.4 = 


y^ Qis-R.) 0isn) 


(4.9) 


{E (E 

{■5^} {«7j} 

The variance of A has been computed using the standard techniques of autocorrelation 
analysis [^, Since the autocorrelation function has a very long tail (due to the 
fact that r^xp ~ N/f 3> Tint,x where / is the acceptance fraction of the algorithm 
and X a generic 5^/060/observable), the self-consistent windowing method proposed in 
28| does not work. Indeed even using a large window of 40ri„t,^, the autocorrelation 
time is largely underestimated. To get a reliable estimate of Tint we have used the 
recipe proposed in the Appendix of . We compute Tint x by 


M 

rint,x{M) = - + ypx(t) 

^ t=l 

rint,x = T\nt,x{M) + px{M)M log 


' N ' 


(4.10) 

(4.11) 


where px{t) is the normalized autocorrelation function, N the length of the walk, / 
the acceptance fraction; M is determined self-consistently and is the smallest integer 


such that 2Qt' int,x{M) > M. We have checked the ad hoc dehnition ([4.11|) in the case 
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standard pivot 

improved pivot 

geometry 

observable 

'^int 

'^int 

1 

Oi 

6.92 ±0.23 

1.4869± 0.0030 

Pn 

282 ± 60 

5.381 ±0.013 

2 

Pn 

280 ± 60 

5.026 ±0.011 

3 

Oi 

7.73 ±0.27 

1.4946 ±0.0030 

Pn 

307 ±68 

5.081 ±0.012 

4 

Pn 

311 ±69 

4.2485 ±0.0083 

5 

Oi 

8.58 ±0.32 

1.4293 ±0.0029 

Pn 

312 ± 70 

4.628 ±0.010 

6 

Pn 

261 ± 53 

2.7910 ±0.0048 


Table 1: Autocorrelation times for the pivot algorithm and our improvement for the 
observables Oi (in the geometries where it does not vanish) and pti. Here N = 1000. 


of ordinary random walks for which exact results are available |^]: the error in the 
estimation of the tail turns out to be at most 10%. We have thus set the error bars 
on the autocorrelation times adding to the error in the determination of r' int,x{M), 
one tenth of the last term in (^4.11|) . 


4.2 The algorithm 

We have simulated SAWs of hxed length N in three dimensions without any excluded 
region TZ. The walk is given by A^ + 1 lattice points {a;,}, i = 0,..., A^, and always 
starts from the origin, i.e. (Uq = 0. Since all our observables are completely symmetric 
under the cubic group, it is not restrictive to fix the first step; we have chosen uji = 
( 1 , 0 , 0 ). 

The simulation used the pivot algorithm |30| , In the standard implementa¬ 
tion a site ujk in the walk and an element g of the lattice symmetry group are chosen 
randomly, with uniform probability. The proposed configuration, obtained from the 
actual one by applying g to the part of the walk subsequent to cufc, is accepted when¬ 
ever self-avoiding. This algorithm is extraordinarily efficient for the study of global 
observables like, for example, the end-point position; indeed, the integrated autocor¬ 
relation time for these observables grows with the number N of steps in the walk like 
N^, where p ~ 0.11 in d = 3, while the CPU-time to produce an independent walk 
scales like N which is the optimal situation. 

However in our case not all the observables are of global character: an example 
is the probability of intersection p-ji. Indeed, as can be seen from table |^, already at 
N = 1000, Tint,pT^ ~ 300. It is easy to understand the origin of these autocorrelation 
times: indeed suppose the walk intersects IZ and let a be the smallest integer such 
that uja € TZ. Then all subsequent walks in the simulation will also intersect TZ at 
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least until a successful pivot move at point i with i < a is performed. The probability 
of such a move is af /N where / is the acceptance fraction. Now the problem is that 
the typical a is very small: if one considers the set of walks that intersect at least 
one of the six needles of the geometry 6, the mean value of the smallest a such that 
cUq, G 7^ is ~ 13, with very small iV-dependent corrections. Thus a rough estimate 
of the autocorrelation times should be iV/(13/). For N = 1000 we have / ~ 0.45 so 
that we expect r ~ 200 which is indeed the order of magnitude we hnd. Moreover 
we expect r to increase as as it should for local observables. Indeed we hnd 
that r for p-ji increases as with p ^ 0.9 except in the case of geometry 6 where 
we observe p ~ 0.6, which could however well be an effective exponent in the region 
500 < iV < 32000. 

The other observables, like Oi, are instead of more global character and indeed 
the autocorrelation times are much smalleif]. They also increase more slowly with iV; 
for Oi we hnd Tint,Oi ~ 

To eliminate these long autocorrelation times one must increase the frequency 
of the moves with the hrst points of the walk as pivots. We have thus modihed 
the algorithm in the following way: one iteration consists now of the following three 
moves: 

• one move with pivot point in uji. 

• one move with pivot point uJk with k uniformly chosen in the interval 1 < fc < 13; 

• one move with pivot point with k uniformly chosen in the interval 13 < fc < 
N -1. 

With this improvement, the autocorrelation times, are greatly reduced. Table 
contains a comparison between autocorrelation time for the two algorithms for walks 
of length Ntot = 1000. We observed sensible reductions of autocorrelation times for 
each observable, but the most impressive results are obtained for the observable 
which denotes the fraction of walks that intersect the considered excluded set. To 
compare CPU-times one should notice that one iteration of the improved algorithm 
takes three times the CPU-time of an iteration of the standard algorithm. Thus in 
practice use of the improved algorithm allowed ns to gain a factor of 15 — 20 on 
p-ji (although no improvement on the exponent p) and a factor 1.5 — 2 on global 
observables like Oi. 


5 The results 


We have studied SAWs with length 500 < N < 32000. The number of iterations for 
each value of N is reported in table The total simulation required 2500 hours of 


^When we speak of autocorrelation time for Oi we really intend the Tint,A as defined in (4.9) as 
this is the observable which controls the error on Oi. 
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500 

1000 

2000 

4000 

8000 

16000 

32000 

^iter. 

3 • 10^ 

5 • 10^ 

6 • 10^ 

o 

r—1 

11- 

1 • 10^ 

5 • 10^ 

2 • 10^ 


Table 2: Number of iterations. 


CPU on a AlphaStation 600 Mod 5/266. We have computed A studying the quantities 
we have previously discussed, whose leading behaviour for A —>• cx) is 

{Ok)N = Bn{k)N- (5.1) 


where cr = fcz/ — A, /c = l,2,4. We have performed standard power-law hts neglecting 
the next subleading terms. This introduces a systematic error which in our case could 
be particularly serious due to small value of A. Indeed one expects corrections to ( |5.1| ) 
of the form which decay very slowly and could thus give sizable corrections 

even at the relatively large values of N we use. To get an idea of the systematic error 
we have repeated the hts using each time only those values of N satisfying N > Ncut- 
in this way one gets different estimates which should converge to A as Ncut goes to 
inhnity. If the different estimates are essentially independent of Ncut (within error 
bars) one can reasonably trust the estimate of A, otherwise one gets an idea of the 
size of the systematic error. 

We show the results of our hts for each geometry in Tables In the hrst column 
one can hnd the diherent values of Ncut-, in the second one the raw Monte Carlo data, 
in the next one the estimated value of the exponent and hnally the conhdence level 
of the ht. In the last column we report the “classical prediction” for the exponent a, 
obtained using (|2.8|) . In the computation of the expectation values we used the 
estimate for v given by p7[ 

= 0.5877 ±0.0006 (5.2) 


Let us now comment on the results. Let us consider hrst the observable Oi. Since 
the estimates for Ncut > 4000 are essentially constant within error bars one could 
estimate 

r 0.396(3) for geometry 1 

(T = < 0.388(3) for geometry 3 (5.3) 

[ 0.379(3) for geometry 5 

These estimates have very small statistical error bars. However the fact that the 
estimates do not agree within the stated errors is a clear indication that there is a 
much larger error due to the neglected corrections to scaling. Indeed a closer look to 
the data shows that in all cases there is still an upward trend, although not statistically 
signihcant as the change of a is smaller than the error bars. It is thus more cautious 
to interpret the results for a as lower bounds on the true value. Thus we estimate 
a = 0.39 ^ 0.01 and thus A < 0.20 ± 0.01. 

Let us now consider the observable 02 - This quantity shows much larger correc¬ 
tions to scaling compared to the previous one: indeed in no case we can identify a 
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geometry 1 


cut 

(0)±W 

a ±Sa 

CL 


Oi 

500 

-1.5259 ±0.0014 

0.40044 ± 0.00044 

9 • 10“^% 

u-A = 

l-u = 

0.4123 ±0.0006 

1000 

-2.0229 ± 0.0018 

0.39831 ± 0.00060 

0.114% 

2000 

-2.6751 ± 0.0028 

0.39584 ± 0.00091 

14.118% 

4000 

-3.5322 ± 0.0044 

0.3932 ± 0.0016 

62.627% 

8000 

-4.6324 ± 0.0060 

0.3955 ± 0.0029 

92.618% 

16000 

-6.092 ±0.015 

0.3962 ± 0.0080 


32000 

-8.018 ± 0.040 



O 2 

500 

-20.457 ± 0.047 

0.97813 ± 0.00098 

29 • 10“^% 

■ 

1000 

-39.944 ± 0.089 

0.9818 ±0.0013 

3.622% 

2000 

-78.25 ±0.20 

0.9859 ± 0.0020 

53.836% 

4000 

-155.00 ±0.43 

0.9866 ± 0.0031 

35.367% 

8000 

-306.42 ± 0.85 

0.9908 ± 0.0054 

26.964% 

16000 

-606.0 ± 3.0 

1.005 ±0.014 


32000 

-1217. ± 10. 



O 4 

500 


2.277 ± 0.035 

15.228% 

4i/-A = 
l + 2u = 

2.1754 ±0.0012 

1000 


2.239 ± 0.044 

20.112% 

2000 


2.182 ±0.063 

22.338% 

4000 


2.056 ± 0.096 

51.065% 

8000 


2.10 ±0.20 

25.931% 

16000 


1.43 ±0.62 


32000 





Table 3: One excluded needle. Here CL denotes the confidence level of the fit. 


geometry 2 


cut 

(0)±,5(0) 

a ±Sa 

CL 


O2 

500 

-45.97 ±0.11 

0.9622 ± 0.0010 

22 • 10“^ % 

;—1 

II 

<1 

1 

CN 

1000 

-88.71 ±0.20 

0.9664 ± 0.0014 

1.206% 

2000 

-171.86 ±0.46 

0.9712 ±0.0021 

32.894% 

4000 

-336.7 ± 1.0 

0.9732 ± 0.0034 

23.052% 

8000 

-658.2 ±2.0 

0.9801 ± 0.0060 

33.624% 

16000 

-1292.7 ±6.9 

0.994 ±0.016 


32000 

-2575. ± 25. 



O4 

500 


2.206 ±0.015 

11.869% 

4jy - A = 
l±2z/ = 

2.1754 ±0.0012 

1000 


2.190 ±0.019 

14.368% 

2000 


2.162 ±0.028 

15.761% 

4000 


2.099 ± 0.045 

39.275% 

8000 

(-1374 ±54) • 10^ 

2.122 ±0.085 

18.371% 

16000 


1.80 ±0.26 


32000 





Table 4: Two opposite needles. Here CL denotes the confidence level of the fit. 
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cut 

(0)±5(0) 

500 

-1.6573 ±0.0015 

1000 

-2.1859 ±0.0020 

2000 

-2.8761 ± 0.0030 

4000 

-3.7761 ± 0.0048 

8000 

-4.9246 ± 0.0065 

16000 

-6.436 ± 0.016 

32000 

-8.435 ± 0.042 


geometry 3 


a ± 6a 


0.39267 ± 0.00044 


0.39038 ± 0.00060 


0.38767 ± 0.00092 


0.3850 ± 0.0016 


0.3873 ± 0.0029 


0.3901 ± 0.0081 


CL 


4 • 10-« % 


0.037% 


13.770% 


59.748% 


70.556% 


u-A 
l-u = 

0.4123 ±0.0006 


500 

-10.696 ±0.025 

1000 

-20.740 ± 0.047 

2000 

-40.42 ±0.11 

4000 

-79.86 ± 0.24 

8000 

-157.19 ±0.45 

16000 

-309.9 ± 1.6 

32000 

-619.8 ±5.6 

500 

-0.1650 ±0.0081 

1000 

-0.276 ±0.016 

2000 

-0.572 ± 0.035 

4000 

-0.819 ± 0.079 

8000 

-1.07 ±0.18 

16000 

-2.02 ± 0.67 

32000 

-2.2 ±2.0 

500 


1000 


2000 


4000 

(-214 ± 12) • 10^ 

8000 

(-871 ± 54) • 10^ 

16000 


32000 



0.9717 ±0.0010 


0.9761 ±0.0013 


0.9809 ± 0.0020 


0.9811 ± 0.0034 


0.9859 ± 0.0057 


1.000 ±0.015 


5 • 10 “^% 

2.063% 

55.682% 2u-A = l 
35.520% 

31.462% 


0.772 ± 0.036 


0.730 ± 0.057 


0.512 ±0.092 


0.53 ±0.18 


0.77 ± 0.43 


0.1 ± 1.4 


5.684% 


4.327% 


89.526% 


74.224% 


63.693% 


2z/ - A = 1 


2.244 ± 0.025 


2.211 ±0.032 


2.171 ±0.046 


2.071 ± 0.072 


2.10 ±0.14 


1.56 ±0.44 


10.182% 


16.821% 


17.299% 


40.820% 


18.912% 


4i/-A = 
l + 2u = 

2.1754 ±0.0012 


Table 5: Two needles along two different axes. Here CL denotes the conffdence level 
of the fit. 
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geometry 4 


cut 

{0)±5{0) 

a L 5a 

CL 


O 2 

500 

-25.876 ±0.071 

0.9439 ± 0.0012 


2n - A = 1 

1000 

-49.00 ±0.13 

0.9501 ± 0.0016 

0 . 747 % 

2000 

-93.52 ±0.29 

0.9561 ± 0.0024 

34.436% 

4000 

-182.09 ±0.63 

0.9552 ± 0.0040 

19.884% 

8000 

-351.1 ± 1.2 

0.9643 ± 0.0072 

35.181% 

16000 

-682.1 ±4.1 

0.981 ± 0.020 


32000 

-1347. ± 17. 



O 4 

500 

(-518 ± 13) • 10^ 

2.1760 ±0.0099 

3.475% 

4n - A 
l±2z/ = 

2.1754 ±0.0012 

1000 

(-2425 ± 55) • 10^ 

2.161 ±0.013 

6.562% 

2000 

(-1104 ±27) • 10^ 

2.142 ±0.019 

7.781% 

4000 

(-521 ± 14) • 10^ 

2.094 ± 0.031 

24.042% 

8000 

(-2181 ±61) • 10^ 

2.1083 ±0.057 

9.664% 

16000 

(-991 ±42) • 10^ 

1.86 ±0.16 


32000 

(-359 ± 37) • lO*’ 




Table 6; Four needles in a plane. Here CL denotes the confidence level of the fit. 


geometry 5 


cut 

(0)±,5(0) 

a L5a 

CL 


Oi 

500 

-1.8113 ±0.0017 

0.38443 ± 0.00045 

13 • 10-^*^% 

I/-A = 

l-u = 

0.4123 ± 0.0006 

1000 

-2.3764 ± 0.0022 

0.38193 ± 0.00061 

0.012% 

2000 

-3.1097 ±0.0034 

0.37900 ± 0.00093 

12.449% 

4000 

-4.0587 ± 0.0053 

0.3762 ± 0.0016 

55.622% 

8000 

-5.2613 ± 0.0073 

0.3786 ± 0.0030 

58.046% 

16000 

-6.833 ±0.016 

0.3828 ± 0.0082 


32000 

-8.909 ± 0.046 



O 2 

500 

-0.1892 ±0.0088 

0.763 ± 0.034 

3.242% 

2n - A = 1 

1000 

-0.311 ±0.017 

0.732 ± 0.052 

2.100% 

2000 

-0.645 ± 0.038 

0.518 ±0.084 

81.853% 

4000 

-0.952 ± 0.086 

0.50 ±0.17 

63.392% 

8000 

-1.20 ±0.17 

0.79 ± 0.39 

61.993% 

16000 

-2.27 ±0.66 

0.2 ± 1.3 


32000 

-2.5 ±2.2 



O 4 

500 

(-282 ± 12) • 10^ 

2.211 ±0.016 

5.575% 

4n-A = 

1 ±2n = 

2.1754 ±0.0012 

1000 

(-1393 ±53) • 10^ 

2.188 ±0.021 

10.689% 

2000 

(-649 ± 26) • 10^ 

2.159 ±0.030 

11.283% 

4000 

(-321 ± 13) • 10“^ 

2.088 ± 0.048 

30.950% 

8000 

(-1334 ±54) • 10^ 

2.103 ±0.089 

12.901% 

16000 

(-618 ±42) • 10^ 

1.73 ±0.26 


32000 

(-205 ± 35) • lO'^ 




Table 7: Three needles along three different axes. Here CL denotes the confidence 
level of the fit. 
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geometry 6 


cut 

(0)±W 

a L5a 

CL 


O 4 

500 

(-787 ± 16) • 10^ 

2.1540 ±0.0086 

4.443% 

Av- A = 
l + 2v = 

2.1754 ±0.0012 

1000 

(-3606 ± 67) • 10^ 

2.141 ±0.011 

7.879% 

2000 

(-1611 ±35) • 10'^ 

2.123 ±0.017 

9.185% 

4000 

(-746 ± 17) • 10^ 

2.078 ± 0.028 

33.022% 

8000 

(-3093 ± 76) • 10^ 

2.097 ± 0.053 

15.301% 

16000 

(-1382 ±56) • 10^ 

1.89 ±0.15 


32000 

(-513 ± 50) • lO*^ 




Table 8 : Six needles along the two directions of all the three axes. Here CL denotes 
the conhdence level of the £t. 


region of N^ut where the estimates are constant. In all cases the estimates of a are 
clearly increasing with Ncut- Using the results with Ncut = 8000 we have 


r 0.991(5) 
^ I 0.980(6) 
^ ~ I 0.986(6) 
[ 0.964(7) 


for geometry 1 
for geometry 2 
for geometry 3 
for geometry 4 


(5.4) 


We conclude thus that a > 0.99 ± 0.01 so that A < 0.19 ± 0.01. 

Consider now the observable O 4 . It has much larger errors than the two previous 
ones. Here we can only give a very rough estimate a ~ 2.10(10) so that A 0.25(10). 

Finally let us discuss the results for O 2 . In this case the estimates of a are much 
lower than expected and indeed for both geometries 3 and 5 they seem to indicate, 
although without much conhdence, a ~ 0.5±0.2. It thus appears that this observable 
does not couple to the leading operator breaking the rotational invariance. This fact 
can be proved rigorously for the ordinary random walk: in Appendix A we show that 
for geometries 3 and 5 we have indeed 


{02)n — 0 


(5.5) 


for all values of N. For the SAW {02 )n 7 ^ 0; however our data show that for this 
observable, in formula (p.7|) , not only A but also B-ji vanishes. Thus the study of 
{02 )n allows to compute a new subleading exponent A' > A; from a ~ 0.5 we would 
get A' ~ 0.7. Notice however that the error bars are too big to really trust this 
estimate. 

Let us now discuss the behaviour of the intersection probabilities p-ji. As dis¬ 
cussed before, pure dimensional arguments suggest that the probability that a SAW 
intersects any one-dimensional set is vanishing in a three-dimensional space. This 
argument works of course in the continuum limit. For SAWs on the lattice this 
statement translates in the fact that given a SAW of length A, and a 
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geometry 


cut 

Pn ± Spiz 

500 

0.182995 ±0.000067 

1000 

0.194002 ±0.000070 

2000 

0.203434 ±0.000089 

4000 

0.21173 ± 0.00012 

8000 

0.21861 ± 0.00014 

16000 

0.22461 ± 0.00027 

32000 

0.23041 ± 0.00058 


A±6A CL 


0.2150 ±0.0064 74.926% 


0.215 ±0.012 58.808% 


0.223 ±0.026 40.621% 


0.16 ± 0.13 56.695% 


500 

0.34624 ± 0.00012 

1000 

0.36415 ± 0.00012 

2000 

0.37928 ± 0.00015 

4000 

0.39235 ± 0.00019 

8000 

0.40309 ± 0.00022 

16000 

0.41235 ±0.00044 

32000 

0.42116 ± 0.00094 


0.2357 ±0.0067 I 83.318% 


0.234 ±0.013 69.572% 


0.237 ±0.027 48.991% 


0.176 ±0.052 61.677% 


500 

0.33456 ± 0.00011 

1000 

0.35209 ± 0.00011 

2000 

0.36696 ± 0.00014 

4000 

0.37993 ± 0.00018 

8000 

0.39061 ± 0.00021 

16000 

0.39983 ±0.00041 

32000 

0.40874 ± 0.00088 


0.2276 ±0.0063 I 74.457% 


0.226 ±0.012 58.934% 


0.232 ±0.026 39.799% 


0.168 ±0.086 54.075% 


500 

0.58302 ± 0.00015 

1000 

0.60509 ± 0.00015 

2000 

0.62332 ± 0.00018 

4000 

0.63875 ± 0.00023 

8000 

0.65123 ± 0.00025 

16000 

0.66182 ±0.00047 

32000 

0.6718 ±0.0010 

500 

0.46027 ± 0.00013 

1000 

0.48113 ± 0.00013 

2000 

0.49865 ± 0.00016 

4000 

0.51380 ± 0.00021 

8000 

0.52619 ± 0.00024 

16000 

0.53683 ±0.00046 

32000 

0.5470 ±0.0010 

500 

0.74526 ± 0.00016 

1000 

0.76491 ± 0.00014 

2000 

0.78069 ± 0.00017 

4000 

0.79389 ± 0.00019 

8000 

0.80432 ± 0.00020 

16000 

0.81312 ± 0.00035 

32000 

0.82126 ± 0.00067 


0.2637 ±0.0067 I 75.629% 


0.259 ±0.013 64.062% 


0.258 ±0.027 43.116% 


0.195 ±0.064 55.552% 


0.2397 ±0.0063 I 72.307% 


0.236 ±0.012 58.015% 


0.242 ±0.025 38.594% 


0.177 ±0.072 53.900% 


0.2897 ±0.0068 I 49.055% 


0.280 ±0.012 49.324% 


0.276 ±0.024 30.554% 


0.206 ±0.052 57.904% 


Table 9: Probabilities of intersection pji for the various geometries. Here CL denotes 
the conhdence level of the £t used to determine A. 
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one-dimensional set S such that ujq is at a distance of order from S, then the 
probability that the walk intersects S goes to zero for N ^ oo a.s N~^. In our case 
Pti is however the probability that the SAW intersects a one-dimensional set TZ which 
is at a hnite hxed distance from the origin of the walk. In this case we expect the 
intersection probability to tend to a constant as N ^ oo and thus a behaviour of the 
form 

bn 


PniN) =pnioo) 


Ata 


(5.6) 


The results for pn for the various geometries are reported in Table ^ together with the 
estimates of A from a three-parameter £t of the form (|5.6|) . The estimates indicate 
A = 0.23(3) except for the geometry 6 in which case one would derive a much higher 
value of A, A 0.28. These discrepancies should not however be taken seriously: 
the fit ( b.6|) is very unstable in presence of additional subleading corrections. To 
understand the size of the systematic error one should expect from fits of the form 
( b.6|) , we have performed the following test: we have considered Oi and we have 
analyzed the data as 


< Oi> 


N 




= a+ hN 


-A 


(5.7) 


where we have used v = 0.5877. In all cases we have obtained estimates A 0.16 — 
0.23 and moreover we have found a barely compatible with zero within error bars (for 
instance for the geometry 1, Ncut = 4000, we have a = (24 ± 23) ■ lO”"^). Clearly the 
additional corrections play still an important role. It is however reassuring that the 
value of A is in essential agreement with what we expect. 


6 Conclusions 

In this paper we have studied the role played by one-dimensional vacancies in the 
critical behaviour of the three-dimensional SAW. As already pointed out in [|l], a new 
critical exponent arises. We have carefully checked that the exponent depends only 
on the dimensionality of the vacancies by verifying its independence from the shape 
of the excluded region: in particular it does not depend on the discrete symmetry 
which the lattice has after the introduction of the excluded region. 

We have given a geometrical argument to interpret the new exponent and we have 
thus derived a classical prediction for it. Of course we expect renormalization effects 
to change the classical formula: we will thus write 

A = 2v — 1 + pn^ 0.175 -I- pn (6.1) 

The quantity pn is an anomalous dimension which is expected to be small. Our 
numerical data give 

A = 0.18 ±0.02 (6.2) 

so that \pn\ ^ 0.02. The classical prediction, obtained setting pTi = 0 is thus a very 
good approximation. 

Let us notice that our estimate of A is somewhat lower than the estimate of [|l|, 
A 0.24 and of [^, A 0.217 ± 0.013. The origin of these discrepancies is probably 
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in the neglected additional corrections to scaling: the results of are obtained from 
an exact enumeration and thus probe only very short walks, while the estimate of [0 
comes from walks which have mainly N ^ 1000 — 4000. Here we use longer SAWs and 
a much higher statistics and thus we can do a much more detailed study of the role 
of the next subleading terms. We thus hope to have a better control of the additional 
corrections although it is conceivable that also our present estimate is systematically 
higher than the true result. We hope that our error bar, which we believe is very 
conservative, takes correctly into account these systematic effects. 

Let us hnally remark that all the arguments we have given are of extremely general 
nature and can thus be used for other systems and geometries. 


A Critical behaviour of random walks in presence 
of dji dimensional vacancies 

Consider a d-dimensional lattice H and a region TZ. Then let cn{x) be the number of 
ordinary random walks, starting from the origin and ending in x, that never intersect 
TZ except for f = 0. In this appendix we will derive a general integral equation for 
the generating function which can be solved exactly when 7?. is a dT^-dimensional 
hyperplane. In this way we will be able to check the computations of Section 2. 

Let us start from the recursion relations 

Co (a:) = 6 (x, 0) , 

cn {x) = [civ-i + cjv-i (x - e*)] ■ 1 - ^ d(x, a) 

i=l \ 

where e* is the unit vector in the ^direction. Let us now introduce the generating 
function 

CXD 

G{x) = I3 ^cn{x) (A.2) 

Af=0 

and its Fourier transform G{q). It is then a simple matter to obtain the following 
equation 

G{q) = 1j3G{q) {2d — 

“'5 /' 7VV - *") L (a.3) 

J-TT (27r) 

where qi = 2 sin (f-)- Dehne now the free propagator which coincides (apart from a 
factor (3) with the generating function for unconstrained random walks 

1 - 2/3 Y cos Qi 


for V > 1 

(A.l) 
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where ttiq = Then ( A.3 ) can be rewritten as 


^ 1 PIT ^ ^ 

G(q) = -D(q) -D(q) G(k)(2d -e)Y. (A.5) 

P J-1T t^TTJ 

This equation is completely general and applies to any excluded region 71. 

Let us now restrict ourselves to the case in which 71 is the dT^-dimensional hyper¬ 
plane given by the equations = ... = Xd = 0. If we now multiply (|A.5|) by 

{2d — (f) J2r)e:'ji and integrate over q, using the identity 2d — (f‘ = ^ — D{q)~^^ 

we obtain 



d^q 

{2tiY 

1 


G{q){2d-f) = 


i-p 


d'^q 

l-qr {2ny 


-1 


D{q) E 


p{q-k')-a 


a&Tl 


Inserting back in ( |A.5|) we get hnally 

G{q) = 


D(q) 


/-I («i. ■ ■ ■. qd^Xd^+u ■■■,kd) 


(A.6) 


(A.7) 


where V-ji = d — dji . It is easy to check that this solution has the correct properties. 
Indeed one can verify immediately that G{x) = 5^^ whenever x E 7Z. Notice hnally 
that for d-ji = 0 the solution does not become the standard solution for unconstrained 
random walks, but gives the generating function for random walks with the excluded 
origin. 

Let us now consider the critical limit mo —>■ 0 and let us dehne 




d^^k 

{2'k)'^'^ 


D{qir--Adn^ 


k 


dqz + l^ ■ 


■■,kd) = 



cP^k 1 
{2t^P^ ml + P 


(A.8) 


where ml = m^ -|-g| and q\\ = (gi,..., qd.j^, 0,..., 0). The introduction of the excluded 
region will thus be relevant or irrelevant depending on the behaviour of /(mg) for 
mo —>■ 0. If the limit exists the perturbation is irrelevant, while if the integral diverges 
the perturbation is relevant or marginal. This last case corresponds to 77-ji = 1, 2. 
For 77-ji = 1, we have 



(A.9) 


and the critical behaviour is drastically changed since 


G{q) 


^ml + gf^mg + + 4 

ml + (p 


(A.IO) 
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In the critical limit we have 


G{q) 


2^ml + qf 

ml + 


From this expression it is easy to verify that, for N ^ oo, 

cn = J2(^n{x) ^ -^={2dfN-^/^ 


(A.ll) 


(A.12) 


so that 7 = 1/2 which differs from the value of 7 for random walks in free space, 
7 = 1. Analogously we have 


- 21 N 

Xi)n ~ -J , 

(A. 13) 

2AA 

- -J- • 

(A. 14) 


The exponent u is not changed, but the amplitudes are, as expected, dependent on 
the direction and different from the value they assume in free space, {xDn = {xI)n = 
N/d. 

If Vti = 2, which, as we shall see, corresponds to a marginal operator, we have 



2 1 
Trrho + 4 


K 


4 

rho + 4 


(A,15) 


where K{z) is an elliptic integral. For ml 0, lTi{ml) diverges logarithmically so 
that we get in the critical limit 


G{q) 


Att 

{ml+ q^)log[{ml + qj)/32] 


(A.16) 


Thus the propagator differs from the unperturbed one only by a logarithmic correc¬ 
tion. From this expression we easily get 


Cn = '^CNix) ~ 

X 


{2d)^ 

logiV 


(A.17) 


Thus in this case we have 7 = 
tional logarithmic correction 
we hnd 


= 1 as for random walks in free space: however an addi- 
appears as expected in the marginal case. Analogously 





N 

logiV 


(A. 18) 


If now Vti > 2 I-ji{ml) has a hnite limit for m^ —0 so that in the critical limit 


G{q) 


Djq) 

/7^(0) 


(A.19) 


23 



Thus G{q) is identical to the generating function of unconstrained random walks 
except for a multiplicative constant which is related to the total number of walks 
ctv = J2x cn{x). Indeed from ( |A.19|) it follows that for large N 


{2d)^ 


(A.20) 


so that pti, the probability that an unconstrained random walk intersects TZ is simply 
1 — l/( 2 d/ 7 ^( 0 )). It is easy to check that this probability tends to zero as d —>• cx) at 
(Iti hxed. Indeed, let us compute the large-expansion of we start from the 

standard representation of in terms of Bessel functions: 

In{ml)= r (A.21) 

Jo 

where Io{t) is the zeroth-order Bessel function of first kind. Expanding Io{t) around 
f = 0 we get hnally 




1 

ml 2'D'ii 

(2?7-i H-h _1_ iA 221 

n^o ni"=o (tjiI + 2V^^fGl+■■■+nT,^) 


so that, for V-ji —>• oo, 


ln{0) 


1 



1 



(A.23) 


It follows that for d —i> cx), ptj ~ {2dTi + l)/(2d). 

Let us now discuss the subleading corrections to ( |A.19|) . We need here the small- 
tuq expansion of J(mg). We start from the well-known asymptotic expansion for large 
t of the Bessel function Io{2t) 






=0 {2ty 




Using (|A.21|) we can rewrite our integral as 


(A.24) 




dte-^^i [e-^*Jo(2t) 




dte-^^o 


+ 




bn(JJn) 

-+n 


=0 (2t)- 


+ 


dt 


(A.25) 


The hrst two integrals have a regular expansion in terms of powers of ml 
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Let us now determine the behaviour in of the generic term appearing in the 
series in the last term. Integrating by parts we get 


' " 9 


bnjT^n) 


^VT^+2n-2 

‘2'D-Ti/2+n 


bn{Vn) di 




/mo 


^Vn/2+r 


g{k,n) 

E 

i=i 


■D'-'n 


(-1 


f g{k,n) 

n 


i=i + n-i^ 

1 


E- 

. 5=0 


-ir 
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s'. 


+ 


=1 'DTi/2 + n-i 


(A.26) 


where, for k integer, 
g{k,n) = 


k + n — 1 if V-ji = 2k] 
k + n + 1 if V-ji = 2k + 1 

Ei(-hio) = - log mg - - Y.T=i = 2 ^; 


E("io) - 1 v^_ - 3y-oo (-i)'=m; 


- < Er=i 


0 (fc+f).fc! 


if P 77 = 2 fc + 1 


(A.27) 

(A.28) 


Here 7 ^ is the Euler-Mascheroni constant, 7 ^ ~ 0.5772156649. The whole integral 
can be represented in terms of mg by 


Inifnl) = E ^nhig” + mQ^ ^C'(mg) ; 


(A.29) 


n=0 


An are suitable constants and C'(mg) is a function of mg which is hnite for mg —0 for 
V-ji odd, diverging logarithmically for T>ti even. The second term in (|A.29| ) represents 
the effect of the excluded region and corresponds to an exponent A 


A — v{Vti — 2) — -^n — 1 


(A.30) 


which agrees with our prediction. 

To conclude this appendix let us prove a result for more general excluded regions 
which we will use in the main text. Let the excluded region be of the form TZ = 
where 77.* is some subset of points of the Ath coordinate axis. If A 7 ^ A • • • ^ im 
2 < n < d, then 

{xi^Xi^.. .Xi^)N = 0 (A.31) 

In particular, for geometries 3 and 5 we have {02)n = 0. 

To prove (|A.31|) consider (|A.5|) , which, as we already said, is valid for general 
excluded regions 77. Using now 


d d 


d 


dqi^ dqi^ dqi 

d d d 


■D{q) 


= 0 


q=0 


dqi, dqi^ dq, 


Qi{k-q)-a _ g 


(A.32) 

(A.33) 
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N 

{x^ + y^ + z^)N 

{x^ + y^ + z^)n 

Qn 

500 

1000 

2000 

4000 

8000 

16000 

32000 

1785.8 ± 1.1 

4051.8 ±2.0 
9167.4 ±4.4 

20752.9 ±9.4 
46995 ± 18 

106166 ± 60 
240356 ± 221 

(28727 ±38) • 10^ 
(14826 ± 16) ■ 10^ 
(76027 ±81) ■ 103 
(38994 ± 39) ■ 10^ 
(20020 ± 17) ■ 10*^ 
(10222 ± 13) ■ 10® 
(5247 ± 11) ■ 10^ 

0.90078 ± 0.00046 

0.90309 ±0.00039 

0.90464 ± 0.00038 

0.90539 ±0.00036 

0.90649 ±0.00032 
0.90692 ±0.00047 
0.90829 ±0.00078 


Table 10: Mean values for SAWs in absence of any excluded region, (a:, ?/, z) are the 
coordinates of the end-point of the walk. 


for A 7 ^ ^2 • ■ • 7 ^ ^n, 2 < n < d, we get 

d d d 


dqi, dqi^ ''' dqi„ 

from which ( |A-31|) immediately follows. 


G{q) 


= 0 


q=0 


(A.34) 


B Large-distance behaviour of the two-point func¬ 
tion 


In this Appendix we will present some results which concern the three-dimensional 
self-avoiding walk with no excluded region and we will use it to discuss, along the 
lines of |2^, the behaviour of the two-point function G(r; (3) in the large-distance 
region |r| > Re{P) where Re{P) is the mean end-to-end distance. Consider now the 
Fourier transform G{p; /3); for /3 —> /3c = 1/p standard scaling theory predicts 


G{p; (3) 
G(0;/3) 


G{q) 


(B.l) 


where q = pRe{(3)/6. An important characteristic of G{q) is the fact that in the 
region < 1, G{q) is essentially a free-held propagator, i.e. it can be parametrized 
as 

G(,) « ^ . (B.2) 

The deviations are small and can be parametrized by a term, i.e. by 


G{q) 


1 

l + q"^ + b2{q^Y 


(B.3) 
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A strong-coupling (exact-enumeration) study set a bound on 62 : 

- 3 - 10 -^< 62<0 . 


(B.4) 


The constant 62 has also been computed [^] in the e-expansion and in the expansion 
in fixed dimension with the result : 62 = ~3 ■ 10 “^. 

Here we want to give a bound on 62 using our Monte Carlo data. A simple 
computation gives 

1 ri'-v'iFi'o/ -u AiA 

(B.5) 


i r(7)r(7 + 4z/) ^ 
^ 2 r(7 + 2z/)2 ^ 


where 


Q = lim Qm = lim 

N^oo N 


+ + z^) 


N 


00 (x^ + y'^ + z‘^)% 


(B. 6 ) 


Our Monte Carlo estimates for are reported in Table |^. It is evident that the 
data show strong corrections to scaling. To determine Q we have thus performed a 
fit of the form 


We find 


Q 

A 

2 


X 


-«+;VA ■ 

(B.7) 

0.9091 ±0.0016 , 

(B. 8 ) 

0.41 ±0.16 , 

(B.9) 

1.77 (4d.o.f.) . 

(B.IO) 


The value of A is in agreement with the estimates of |^. Using for 7 the value p2 
7 = 1.1575(5) we get finally 


62 = -(13± 17) 10 


-4 


(B.ll) 


Our Monte Carlo data confirm the fact that &2 is extremely small although we are 
unable to compute the actual value. 

together with the estimates of 7 


On the other hand we can use (p.4|) and 


and u to obtain an estimate of Q. We get Q = 0.9082(11). 
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